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Abstract 

Matrix pair beamformer (MPB) is a promising blind beamformer which exploits the temporal sig- 
nature of the signal of interest (SOI) to acquire its spatial statistical information. It does not need 
any knowledge of directional information or training sequences. However, the major problem of the 
existing MPBs is that they have serious threshold effects and the thresholds will grow as the interference 
power increases or even approach infinity. In particular, this issue prevails in scenarios with structured 
interference, such as, periodically repeated white noise, tones, or MAIs in multipath channels. In this 
paper, we will first present the principles for designing the projection space of the MPB which are closely 
correlated with the ability of suppressing structured interference and system finite sample performance. 
Then a multiple-interference-channel based matrix pair beamformer (MIC-MPB) for CDMA systems is 
developed according to the principles. In order to adapt to dynamic channels, an adaptive algorithm for 
the beamformer is also proposed. Theoretical analysis and simulation results show that the proposed 
beamformer has a small and bounded threshold when the interference power increases. Performance 
comparisons of the MIC-MPB and the existing MPBs in various scenarios via a number of numerical 
examples are also presented. 
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Adaptive arrays, code division multiple access (CDMA), matrix pair beamformer, structured inter- 
ference. 

I. Introduction 

Adaptive beamforming is a promising technique to spatially suppress interference, and can be used 
in dense interference environments, such as, direct sequence code division multiple access (DS-CDMA) 
systems. Adaptive beamforming techniques often make use of a known training sequence or the direction- 
of-arrival (DOA). However, the time-varying nature of mobile communication requires continuous DOA 
tracking or pilot signals in these methods, which increases the complexity and bandwidth requirement. 
In addition, steering vector errors will cause performance loss in DOA-based beamformers as well 12, 
0. 

To overcome these problems, many blind adaptive beamforming algorithms have been extensively 
studied. The constant modulus algorithm (CMA) is a class of gradient-based algorithm that works on 
the premise that the existence of an interference causes fluctuation in the amplitude of the array output, 
which otherwise has a constant modulus 0-116]]. But for the possible presence of constant modulus (CM) 
interfering signals (e.g. MAI, BPSK jamming, etc.) and the requirement for power control, the blind 
algorithm based on CM property is less feasible for DS-CDMA systems lfl2ll . Another class of blind 
algorithms exploit the temporal signature of the signal of interest (SOI) to acquire its spatial statistical 
information, which also only requires the spreading code and timings of the desired user f71l- lfT2l as 
the CMA methods |4]-|6]. In iTTTl — ifTDll . the eigenstructures of the pre- and post-correlation (PAPC) array 
covariance matrices are used to derive the beamformer, and various kinds of low complexity iteration 
algorithms are developed. The Maximin algorithm proposed in ifTTI . lfl2l uses a filter pair (FP) to separate 
the SOI and the interference, and update the weight vector by steepest decent method. 

As indicated in our recent work Ifl3ll , Ifl4l . these approaches share the same processing structure, 
i.e., two projections to construct two estimated matrices followed by a generalized eigen-decomposition 
of the matrix pair, and hence are referred to as matrix pair beamformer (MPB). We also find the key 
assumption that the two matrices share the same interference statistics is not valid in many cases, which 
will cause so-called matrix mismatch fl3l . lfl4ll . Due to matrix mismatch, the MPB always suffers from 
a threshold effect. When the input signal-to-noise ratio (SNR) is below the threshold, the performance 
of the beamformer will degrade rapidly, and the main beam will point to the direction of interferers. In 
some cases, the threshold SNR is infinity and the MPB fails forever. Furthermore, the existing MPB is 
vulnerable to structured interference in many cases, such as periodically repeated white noise, tones, and 



October 4, 2010 



DRAFT 



3 



MAIs in multipath channels. As a result, the threshold will grow as the interference power increases. In 
order to make the beamformer work, the power of the SOI should also increase to compete with that 
of the interference. This property means the MPB cannot function under this condition. Therefore, it is 
important to design an MPB with ability of suppressing structured interference. 

Finite sample effect is another important factor having an impact on the performance of a beamformer. 
Since insufficient sample-support may cause a considerable mismatch between true and sample covariance 
matrices in practical implementations, the calculated noise eigenvalues will be a significant spread around 
the correct values [1]. As a result, how much independent noise samples obtained can determine the 
performance of a beamformer. Robust design of a beamformer involving diagonal loading factor IPT51 , 
lTT6l is another approach to cope with this problem, which desensitizes the system by compressing the 
noise eigenvalues of the correlation matrix so that the nulling capability against small interference sources 
is reduced |[T6l . However, how to choose the best loading factor in a real scenario in order to combat 
the finite sample effect is still an open problem. 

Based on the above observations and the analytic results in our recent work, in this paper, we first 
propose several principles for designing the projection space for MPBs. Then a multiple-interference- 
channel based matrix pair beamformer (MIC-MPB) for CDMA systems is developed. The beamformer 
has a small and bounded threshold, i.e., the threshold does not grow when the power of the interference 
increases. Moveover, by exploiting more signal-free interference samples, the approach achieves a less 
perturbed noise subspace and avoids signal cancelation. 

The rest of the paper is organized as follows. Section II presents a general framework of MPB to 
summarize and reinterpret the basic ideas in f7l- |[T2l . followed by reviewing some results concerning 
threshold effects of the existing MPBs. In Section III, we first present the principles for designing the 
projection space based on the results. Then, a multiple-interference-channel based MPB is proposed 
according to the principles. In order to adapt to dynamic channels, Section IV derives an adaptive 
algorithm for the proposed beamformer. Finally, Section V gives a number of computation and simulation 
results that illustrate the good performance of this beamformer, and Section VI concludes the paper. 



II. Problem Formulation 



A. Signal Model 



In a CDMA system with M users, the transmitted baseband signal of the ith user is 

+00 




(1) 



k=— 00 
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where Pt is the transmit power; bi(k) S {+1, —1} is the fcth transmitted symbol by the ith user; a(t) 
is its normalized signaling waveform, supported on [0, T s ]; and T s denotes the symbol interval. a(t) can 
be expressed as 

N-l 

d{t) = Ci(n)^{t - nT c ) (2) 

n=0 

where Ci(n) G {+1,-1} is the spreading code assigned to the ith user; ip(t) is the normalized chip 
waveform with time duration T c ; and N = T s /T c is the processing gain. 

The receiver has an antenna array of L isotropic elements that receives signals from far field. Each 
user signal arrives at the array via different paths. We assume all elements experience identical fading 
for each path. In addition, there are Q jammings received. Then the total received signal after carrier 
demodulation is 

M-1 Di Q 

x(t)= ^^a ij Si(t-Ti j )a(9i j ) + ^2z q (t)aL(e q )+v{t) (3) 

i=0 j=l 9=1 

where a^-, Tjj and a(9ij) are the path gain, delay and array response vector for the jth path of the ith 
user; Di is the number of paths for the ith user; z q (t) and a(9 q ) are the waveform and the array response 
vector for the qth jamming; v(t) is the space-time white noise. For uniform linear array (ULA) with 
interelement spacing d and carrier wavelength A, the Ith component of a(9) is e~ j sm ( d ) , where 9 is 
the DOA and can be 9 q or 9ij. 

After matched filtering and chip-rate sampling, the discrete signal can be written as 

Mn+1)T C 

x(n) = / x(t)ip*(t-nT c )dt 

JnT 

M-1 Di +oo Q 

= 2E h{k)ci{n - riij - kN)a(9 i:j ) + z q (n)a(9 q ) + v(n) (4) 

i=0 j=l k=—oo q=l 

where (•)* denotes conjugate; Py and riij are power and chip delay for the jth path of the ith user, 
respectively. We have omitted aij in (0]) and contained it in P^; z q (n) and v(n) are the discrete counterpart 

of z q (t) and v(i). 

We also assume the propagation delays of multipath signals from a desired user, enumerated as i = 
in (01), can be perfectly estimated as the existing MPBs l71l- |[T2l . and our goal is to recover bo(k) from 
x(n) with fidelity. There are Dq paths for the desired user, and our strategy is to construct beamformer for 
each path to suppress all other signals except the specified path. In fact, the delayed replica of the desired 
signal in the multipath propagation can be treated as MAIs when the relative delay between a certain 
path and the desired one is greater than one chip, since the spreading code is assumed to have good 
cross-correlation and self-correlation property. Then, a two-dimensional rake combiner is employed to 
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combine outputs of the Dq beamformers, and the procedure is similar to |9), ifTOl . Since the main puipose 
of this paper is to address the problem of the threshold effect of the MPB, without loss of generality, the 
first beamformer (corresponding to the first path of the desired user) is used for the following analysis 
for notational convenience. To be more specific, we rewrite © as 

D 

x(n) = VPiSi(n)a(9i) + v(n), (5) 

where D = ZZo 1 A + Q - K L; Si(n) is the discrete sequence of the ith signal with normalized 
power, with so(n) is the SOI, and si(n), S2{n), . . . , sj}(n) are interferers such as other multipath signals 
of the desired user, MAIs by other M — 1 users, and jammers, etc. Pj, a(0j), and 0+ are its power, steering 
vector and DOA, respectively. Specifically, the SOI so(n) is 

+ 0O 

s o{n) = ^ b (k)c (n - kN - n ), (6) 

k=— oo 

where no = noi is the equivalent propagation delay. 

B. The Matrix Pair Beamformer 

The steering vector a(0j) in ([5]) is a spatial signature of the ith signal, which is different from others so 
long as they arrive from different directions. Beamformer is a spatial filter that exploits such difference 
to pass the desired signal so(n) while suppressing si(n) . . . S£>(n) and v(n). A statistically optimum 
beamformer [H generally requires at least, either explicitly or implicitly, the information about the steering 
vector a(#o) and the interference covariance matrix. The latter one may be replaced by the data covariance 
matrix, so the remaining problem is how to acquire sl(9q). To work "blindly", i.e. without explicit 
information of DOA, the methods in iPTl- lfTZl exploit the temporal signature of the desired signal to acquire 
these spatial statistical information. Specifically, it is implemented by two orthogonal projection operations 
and a generalized eigen-decomposition to exploit a "mismatch-match" mechanism in a covariance matrix 
pair. Hence, we refer to them as matrix pair beamformer Ifl3l . |[T4l . With the data segmentation, the array 
outputs corresponding to the kth symbol of the SOI can be expressed in the following matrix form: 

x(kN + n ) • • • *(kN + n + N - 1) 

D 

aoc^+J] VPi^I (k)+V(k) 
i=i 

a c n +A J e|Sf(A : )+V(A;), (7) 



x(fc) = 



Mk) 
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where a« stands for a(0j), (i £ {0,1, ... , D}) and A/ is a matrix whose columns are the steering vectors 
of interferers ai . . . a^; Co is the temporal signature vector of the SOI composed of the spreading code 
and (-) T denotes transpose; V(k) are the matrix form of the noise; Sj(fc) are the matrix form of the 2th 
interferer and Sj(fc) is the matrix whose columns are si(fc) . . . so{k), with 



co 



ai a 2 • • • a D 

co(0) co(l) •• 



Si{k) = Si(kN + n ) 



S/(fc) = 
V(fc)4 



si(fc) s 2 (/c) • 
v(fciV + n ) 



co(JV-l) 

• Si(kN + n + N - 1) 

S£)(fc) 

• v(£?iV + n + iV - 1) 



G 



j^diagjp^Ps,--- ,P D }. 



Then, the feth data block in each antenna is projected onto two subspaces: signal space S and 
interference space X, respectively. S is a one-dimensional space with base vector = c / V~N, and T is 
a specifically designed rj-dimensional space with base vectors h^, . . . , h^ T \ The projection operation 



hi 1} hi 2) • • • h 



(n) 



produces signal snapshot x^(A;) and the interference snapshot Xx(fe). Define Hj = 
and assume H-^Hj = I, where denote conjugate and transpose. Then the projection procedures 
may be written as 



xs(*0 = X(A;)h* 5 



a + ^=A / e|Sf(fc)cS+v 5 (A;) 
V TV 



Xx(fc) = X(A;)HJ 

'P o &o(*0l a c r H I +A / e|Sf(fc)Hj+V x (A;), 



(8) 



(9) 



where v s (k) = V(k)h% and V x (k) = V(A;)H*. 

Assume the SOI is uncorrelated with the interferers, we can derive the covariance matrices of xs(fc) 
and Xj(fc) as 



R 5 = ^{x lS (A;)xf = o-i aoa? + Q s 
R X 4 l^{x x (fc)X^ (A;) | = a^aoaf + Qx, 



(10) 
(11) 
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where 

a 2 So = P c$P s c = NP (12) 

4 = ^<#Pjc , (13) 

Qs and Qj are the covariance matrices of the last two terms in ([8]) and ©, respectively. P5 and Pj 
are the projection matrices of S and X, defined as 

P 5 = hsh§ = ^c c^ (14) 
Px = H x Hf = f;4 r) [hf ) ] ff . (15) 

r=l 

In practice, R5 and Rj are computed by sample averaging (c.f. Section [TVT >. 

In most of the existing approaches, X is one dimensional space (rx = 1)- The pre- and post-correlation 
(PAPC) scheme 1171- lfTOl uses x(n) to calculate Rj, thus it is equivalent to selecting one column of 
ItvxTV as Hx, i.e. 



H 7 



••• 1 ••• 



i T 



(16) 



The Maximin scheme in ifTTj and ITT21 employs a monitor filter to isolate the interference, which can be 
interpreted as 



H x = c 



T 

I gj27r/ M F . . . e j2irf MF (N-l) 



(17) 



where /mf £ (0, 1] is the normalized center frequency of the monitor filter (MF), and denotes the 
Hadamard product. 

Under the maximum signal-to-interference-plus-noise ratio (MSINR) criterion, it is well known that the 
optimal weight vector for the first propagation path of the desired user w opt is the generalized eigenvector 
corresponding to the largest generalized eigenvalue of the matrix pair (Rs,Rx), i.e., 

R^Wopt = A max Rrw op t, (18) 

where A max is the largest generalized eigenvalue. Therefore, the MPB can maximize the output signal- 
to-interference-plus-noise ratio (SINR) when w opt is applied to x.$(k), and the output y (k) is 

Vo(k) = w^x^/c) = y s (k) + Vl (k) + y N (k), (19) 
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where 



w o^t a o 



yj(k) = -^w^AjefSf (fc)cS 



1 

y N (k) = w^ t v 5 (fc). 

Then, the final array output after a two-dimensional rake combiner can be written as O, flOll 

*(*0 = X>i,o(fc) s (20) 

where yj, {k) is the jth output of the beamformer corresponding to the jth propagation path, and the 
typical expression of yj t0 {k) can be referred to ( fl9l ). 

C. Threshold Effects Regarding MPB 

Based on the theoretical analysis in |[T3l . lfl4l . A max has the following property: 

Amax ~ max|7o + 1,71 + l|, (21) 

where 

_ L(N -flSNR 
70 L/3SNR + iV ^ J 

is a monotonically increasing function of SNR, and SNR = o 2 s Ja 2 is the SNR of the SOI after 

despreading (or equivalently, input SNR per symbol). (3 is the normalized power leakage ratio (PLR) in 

interference channel defined as 

P a °k = N < = So^lSo. (23) 
Po <t 2 So r x 

7i + 1 is the the largest generalized eigenvalue of the matrix pair (Qs, Qz), which is co-determined by 
the structure and power of interferers as well as the projection spaces of the MPB. It can be derived that 
7i could be bounded if the following expression is satisfied IT31 . lfl4l 

i ± nv I cs ± n V/, (24) 

where 

l^n\Px} 



Vi = span|S/(l)|. 
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where (•) denotes orthogonal complement space, IZ(-) denotes the range space of a matrix; V/ is the 
space spanned by interference sequences and S/(l) are the waveforms of the interferers in the first period. 
The optimal weight vector w opt can be approximated by the following equation |fT3ll , Ifl4l 

Hi ■ R^ao if 70 > 71 
w opt ~ < (25) 

l j u 2 R^ 1 a ei if 70 < 71, 
where a ei is an appropriate linear combination of the steering vectors of interferers ai,a2, . . . ,&d, and 
fix, [12 are the coefficients. The expression of w opt means that if 70 > 71, the main beam of the MPB 
will point to the DOA of the SOI; if 70 < 71, the main beam of the MPB will point to the DOA of the 
interferers. Furthermore, if /3 7^ 0, the beamformer will form a notch in the direction of the SOI because 
Rx contains parts of the desired signal. 

Our work also shows that the existing MPBs are vulnerable to structured interference, such as peri- 
odically repeated white noise, tones, and MAIs in multipath channels for d24l) can hardly be satisfied in 
some cases of those scenarios. For periodical interference, d24l ) can be rewritten as the following flT"3l , 

mi 

Tz[-p Vl H x ]DTz[l> Vl hs}, (26) 

or equivalently, 

^{h Vi H^Hi} D nfev^hs}. (27) 

where Py 7 is the projection matrix of the subspace Vi, Hv 7 is a base matrix of the subspace Vi = 
T^IHvj j = 7^.|S/(1)|. If d26l ) does not hold, 71 + 1 will grow as the interference power increases. 

From the above discussion, we see that the threshold effects of MPBs rely heavily on base matrix Hj 
for the interference space X. Therefore, in the following section, we will propose appropriate methods 
to handle this effect as well as finite sample performance by designing appropriate base vectors for the 
interference space. 

III. The Multiple Interference Channel based MPB 

In this section, starting from the above results, we first present the principles for designing projection 
space for MPBs. 
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Chip matched 
filter 



l-N Serial _l\ FFT 
to Parallel — |/ 



PN code 



Chip matehcd 
filter 



'-Hg>-* 

nT Y 



l-N Serial l\ 

to Parallel 



PN code 




Chip matched 
Alter 



l-N Serial 
to Parallel I — |/ 




PN code 



Fig. 1. Projection operations to separate the signal channel and the interference channels using FFT base vectors. 



A. Principles for Designing Projection Space for MPBs 

1) Ability of Suppressing Structured Interference: Since an MPB can work properly only if 70 > 71, 
70 should be as large as possible for a given SNR. ([22l> shows that 70 is a monotonically decreasing 
function of (3, so (3 should be designed as small as possible. It can also be found from (fTTT) and (l25l ) 
that, if f3 ^ 0, there will be the sample-correlation terms between the SOI and the interference -plus-noise 
in Rj because of finite sample effects. Even if 70 > 71, the sample-correlation terms will cause the 
main-lobe unstable as well as a signal cancellation effect in the beamformer output (H, iPTTI . Therefore, 
(3 should be designed to be 0. With (T23T ). we can easily derive that 



/3 = & IC5 J 



(28) 



On the other hand, 71 should be as small as possible for given power of interference. (1261 ) means the 
subspace spanned by the columns of Hj projected onto Vi must contain the subspace spanned by 
projected onto V/. Since 7e.|H Vl H^h 5 | C V/, <[26]> always holds so long as ^|h Vi H^Hx| = V/, 
which means the columns of HjHy 7 are linear independent, i.e., 



(29) 



This expression shows the subspace T should be properly designed in order that the subspace Vi = 
7?.|S/(1)| does not contain any vector which is perpendicular to the subspace 7^.|Hj| = X. 

2) Improving Finite Sample Size Performance: If j3 = and Rj does not contain any component of 
the SOI, the beamformer can be considered as an Miminum Variance Distortionless Response (MVDR) 
beamformer when 70 > 71 by d25l ), and the performance of the beamformer is degraded mostly by the 
disturbed noise space [ 1 ] and at least K 2L samples of data are needed to maintain an average loss 
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ratio of better than one-half (less than 3 dB) lfl8l . It can be considered that the number of independent 
noise samples available is the number of the effective samples. We now examine the relationship between 
the number of effective samples and X. From Q, the rth column of Xj(A;) can be written as 



4\ k ) 



where vj r (fc) = V(k)[h^]*. Since all elements of V(/c) are Ltd zero-mean Gaussian random variables, 
it can be easily obtained that 

E {vT,rik)^A k ')} = ^ S rr'S kk 'l, (31) 

i.e., the noise component vj r (/c) of different x^(fe) is mutually independent. As a result, the number 
of the effective samples extracted per data symbol is rx, and the total number of the effective samples 
is K ■ rx with K symbols. This result shows that the dimension rx of subspace X determines the finite 
sample performance of an MPB. 



T 

a c 



+ A x e]sJ(k) 



.(0 



+ V X , r (k), 



1,2, 



,rx 



(30) 



B. The Multiple Interference Channel based MPB 

According to (128T ) and (|29l ), we can select the subspace X as the following equation 

X = S 1 - = span-tco}- 1 . (32) 

Since only vectors in spanjco} can be perpendicular to X, there is no vector in 7^.|S/(1)| which is 
perpendicular to X so long as Co ^ 7^|S/(1)|. This condition can be easily satisfied in most cases in 
a multi-user CDMA system. On the other hand, the dimension rx of the subspace X equals to N — 1 
under this condition, then the effective number of samples obtained per symbol is also N — 1, which is 
the maximum value obtained when (3 = 0. 

Specifically, we select the following vector as the the rth (r = 1, . . . , N — 1) base vector of the subspace 

X, 



>) 

^t.mic 



c W^, 



(33) 



where |w^, W^, . . . , 1 j are the base vectors of the Discrete Fourier Transform (DFT), defined 



as, 



W 



N 



1 



J^ 1 



(34) 



Comparing with the Maximin or PAPC method which has only one vector in interference channel (or 
equivalently, subspace X ), this method has N — 1 base vectors, so it can be called Multiple-Interference- 
Channel Matrix Pair Beamformer (MIC-MPB). 
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If we define an L x N matrix 



C = [c c • • • c } H /VN, 



an L x N matrix 



and an N x N matrix 



X H (k) = [xs(k) x?»(fc) ... x 



(N-l) 



(*)], 



W = [W° N Wlf 
it can be easily obtained from ([8]) and (f3Qb 



W 



N-l-i 
N i 



X„(k) = X(fc) 



(1) 

Z,MIC 



h 



(AT-1) 
X,MIC 



X(fc)0C c 



W* 



(35) 



(1351 ) indicates the projection operations implemented by the base vectors defined in (1331 ) are equivalent 
to the procedure illustrated in Fig. [Q The zero frequency outputs of all DFTs generate xs(fc), and all rth 

(r) 

frequency outputs form x. x (k). Mixing with the spreading code flattens the spectrum of the interference 
and noise, making the power evenly distributed on all frequencies. Furthermore, using the DFT base 
vectors for projection operations can be efficiently implemented by Fast Fourier Transform (FFT). 



IV. Adaptive algorithm 

In this section, we derive a blind adaptive algorithm for the proposed MIC-MPB for each signal path of 
the desired user. In order to adapt to time-varying environment, we use the exponentially weighted sample 
correlation matrices Rs(fc) and Rx(fc) instead of R^ and Rj. Then, the recursive update equation for 
the matrices can be written as 



Rs(k) = fiK s (k - 1) + x 5 (fc)xf (k) 
Rx(fc) =/zR x (fc-l)+R^(fc) 



(36) 
(37) 



where 



N-l 



R£(fc) 



N 



r=l 



T £4 r, (*)[4 r> (*) 



H 



and n is a positive constant less than 1. Since the update term in (|37T ) is not rank one, we cannot apply 
Woodbury equality ||20l , ||2T1 to compute its inverse. To solve this problem, let x^ (k) = -x!>j-\k) / y/N — 1 
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and define 



H 



r=l 



R X (A:; t) = fiH x (k - 1) + R^(fc; i). 



(38) 
(39) 



Then we have B£(k) = R^(fc;iV - 1), R x (fc;N _ 1) = R x (fc + 1;0) = Rr(fc), and R£(fc;i) 



Rj(fc;t — 1) +Xj^(fe)[xj''(A;)]- H '. As a result, the following recursive equation can be obtained, 

H 



■M, 



K x (k; t) = fi(t) ■ R. x (k; t - 1) + xj } (fc) (jfe) 



where /j(t) is defined as 



At(t) 



/i t = 1 

1 2 < f < JV- 1 
We then apply Woodbury equality to (1401) and obtain 



(40) 



(41) 



c(k;t) 



[//(t^PCM-i)*?^*) 



l + [/.(t)]- 1 4* } W P(fc;t- i)x^(fc) 



i3 



(*)/ 



P(Jfe;t)= [/i(t)]~M I-c(M) x^O) "\P(k;t-l 



H 



(42) 



(43) 



when t = N — 1, the value of P(&; i) are assigned to P(&) = Rj 1 (fe) and reinitialization is need as the 
following, 



P(jfe) = P(k;N- 1) 

P(jfe + 1;0) = P(Jfe;JV-l). 



(44) 
(45) 



In summary, ([36]>, dSJ, (|42l . (l43l . (|44l . and (|45) complete the update of Rx(fc) and P(fc) = Rj 1 ^). 
Then we can update the weight vector w by power iterations 12T1 : 



w(A; + l) =P(k)Rs(k) 
The details of the algorithm are shown in Algorithm [TJ 



w(k) 



w 



(46) 



V. Simulation Results and Discussions 

In this section, we provide numerical examples to verify the validity of the proposed MIC-MPB scheme, 
and compare the performance of it with that of the PAPC and the Maximin beamformer. In the simulations, 
we assume the transmitted DPSK signal is spreaded by a distinct 31-chip Gold sequence (N = 31) and 
modulated onto carrier frequency of 1 GHz for each user. The data-symbol and spreading sequences are 
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Algorithm 1 MIC-MPB Beamforming Alogrithm 

Rs(0) = <5I where S is a small positive number 
P(0,0) = P(0) = S- 1 ! 
w(0) = [ 1 ■■• ] T 
for k = 1, 2 ... do 

Rs(fc) = AtRsO - 1) + x s (fc)xf (fc) 
for t = 1,2, ...,N- 1 do 
if t = 1 then 
/x(t) = a* 
else 

= 1 

end if 

= x<- t) (fc)/v / A r ~H" 

i + [m^I-^x^WI^p^;*- i)4 t} (fe) 

P(fc;t) = [At(t)] _1 {l- c(fc;t)[x< ; ' ) (fc)]- ff }p(fc;t- 1) 
if t = JV - 1 then 

P(fc) = P(fc;AT- 1) 

P(fc + 1;0) = P(fe; AT - 1) 
end if 
end for 

w(fe + l) = P(fc)R s (fc)^g) I 
j, (fc) = w H (fc)x s (fc) 

end for 



randomly generated for each simulation trial at the rates of 100 kbps and 3.1 Mbps, respectively. Since 
each signal path of the desired user is processed separately by employing the two-dimensional RAKE 
receiver, without loss of generality, we assume the desired user has one propagation path in the first two 
subsections. In the last subsection, we will discuss performance of the proposed beamformer in a special 
case for RAKE processing, i.e., there are multipaths with identical delay of the desired user. 

A. Ability of Suppressing Structured Interference 

Firstly, we study the ability of suppressing structured interference of the beamformers. Three typical 
scenarios-the received SOI with periodically repeated white noise, tones, and MAIs in multipath channels 
are simulated with some specially selected simulation parameters of the interferers. In all the cases, 
we consider a uniform linear array (ULA) with eight omnidirectional antennas (L = 8) spaced half a 
wavelength apart. In these simulations, we also assume that the SOI always arrives from 0° and the 
power of the interferers are always assumed to be equal in each scenario. 

Fig. |2] shows the largest and second largest generalized eigenvalues of the matrix pair of the MIC- 
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+ The Largest Eigenvalue Simulated, INR = 10 dB 
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Fig. 2. The largest and 2nd largest generalized eigenvalues of the MIC-MPB vs. SNR in five tones case. 

MPB with five tones interferers. The tones are assumed to impinge on the array from the directions 30°, 
-50°, -20°, 19°, and 45° with frequency offsets 100 kHz, -300 kHz, 0, 400 kHz, and -100 kHz, 
respectively, with respect to the canier frequency of 1 GHz of the SOI. The simulated eigenvalues are 
obtained by computing the matrix pair and Rj from generated received array signals then using 
eigen-decomposition operation. In order to avoid finite sample effects, 1 million data symbols (K = 10 6 ) 
are used to estimate the covariance matrix pair. Theoretical 70 + 1 is computed by (l22l and 71 + 1 by 
using eigen-decomposition of the matrix pair (Q^Qz). From this figure, we can observe that when 
SNR < —0.6 dB, 70 + 1 < 71 + 1 and the largest eigenvalue of the matrix pair equals 71 + 1; when 
SNR > —0.6 dB, 70 + 1 linearly increases while 71 + 1 remains a constant, the largest eigenvalue then 
switches to 70 + 1. Therefore, the threshold of the MIC-MPB can be considered as —0.6 dB. Since 71 + 1 
of the beamformer remains the same when the power of the interferers or the interference-to-noise ratio 
(INR) increases, the threshold of the MIC-MPB is small and bounded in this scenario. 

Fig. |3]-Fig. [5] show the normalized output SINRs corresponding to the MIC-MPB, Maximin, and PAPC 
scheme versus input SNR in the three scenarios. The normalized output SINR is defined as the output 
SINR of the MPB normalized by the optimum SINR with no interference, given by 



SINR. 
SINRl 



opt 



'O 



(47) 



where 



SINR, 



■ A 
>o — 
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SINR opt = • ||a || 2 • ||c || 2 = LSNR. 

The simulated normalized output SINRs are obtained by using the above equations with simulated 
received signals, and the theoretical values are computed by an approximated piecewise function G(SNR) 
described in lH"3l . 1141 . In Fig. [3l two periodically repeated white noise arrive at 30° and —40°, respec- 
tively. The periods of the interferers are both equal to the duration of a CDMA symbol T s . In Fig. 
@1 there is one incident MAI signal with three-ray multipath delays of 3 chips, 5 chips, and 4 chips 
from directions 30°, —20°, and —50°, respectively. The simulation parameters in Fig. [5] are the same as 
those in Fig. |2] Some points need to be noted that these simulation parameters are specially designed in 
order to give prominence to the threshold effects the MPBs, because the threshold of the Maximin or 
PAPC is very small (far more less than SNR) and the beamformers can be well-behaved in most cases. 
Since G reflects the limiting performance of a beamformer, K = 10 6 symbols are simulated for each 
SNR under given INRs in every experiment to eliminate finite sample effects. However, deviation in 
simulated values still can be seen in the figures when INR = 30 dB and SNR are below the thresholds 
of the proposed MIC-MPB scheme. This phenomenon can be explained by d25l ), i.e., when SNR is 
below the threshold, the steering vectors of the interferers will dominate and the beamformer can be 
considered as an Miminum Power Distortionless Response (MPDR) beamformer, which will receive the 
interferers. Since larger INR means more interference power contained in Rj, more data samples are 
required for "satisfactory" performance 12, |fl9ll . But for the Maximin or PAPC beamformer, things are 
totally different. This is because both schemes employ one dimensional interference subspace X, which 
make independent interferers correlated after projection operation. As a result, the steering vector of the 
interferers contained in Rj is a compound vector, which is different from a ei . Therefore, they can be 
considered as MVDR beamformers when SNR are below the thresholds, and far more less samples are 
needed to maintain stable system performance. 

From the figures, we can find that the proposed MIC-MPB scheme can achieve the optimum SINR 
regardless of the received power of interference in the three scenarios when SNR > SNRj , which 
means the structured interference have been totally filtered under this condition. But for the Maximin 
beamformer, more input signal power is needed for it to reach the upper plateau when the power of 
the interferers or INRs increase. Meanwhile, its limiting performance decreases when INR grows. This 
is because the Maximin beamformer cannot perfectly eliminate the interferers in these scenarios, which 
can be verified by Fig. |7J the Maximin beamformer does not form deep nulls in the direction of the 
interferers. For the PAPC beamformer, we can find that it completely fails in the scenarios. Furthermore, 
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Fig. 3. Normalized output SINRs corresponding to the MIC-MPB, Maximin, and PAPC vs. SNR in two periodically repeated 
white noise case. 

TABLE I 

Input SNR Thresholds of the beamformers in two periodically repeated white noise case 



Matrix Pair 


Input SNR Thresholds SNR T0 (dB) 


Beamformers 


INR = 10 dB 


INR = 20 dB 


INR = 30 dB 


MIC-MPB 


-0.93 


-0.85 


-0.84 


Maximin 


7.7 


17.5 


27.5 


PAPC 


oo 


CO 


oo 



TABLE II 

Input SNR Thresholds of the beamformers in three-ray multipath MAI case 



Matrix Pair 


Input SNR Thresholds SNR T0 (dB) 


Beamformers 


INR = 10 dB 


INR = 20 dB 


INR = 30 dB 


MIC-MPB 


-9.4 


-9.3 


-9.3 


Maximin 


6.2 


15.8 


25.8 


PAPC 


oo 


00 


oo 
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Fig. 4. Normalized output SINRs corresponding to the MIC-MPB, Maximin, and PAPC vs. SNR in three-ray multipath MAI 

case. 
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Fig. 5. Normalized output SINRs corresponding to the MIC-MPB, Maximin, and PAPC vs. SNR in five tones case. 



its normalized output SINR decreases to zero in the order of 0(SNR -2 ) when SNR goes to infinity. 

Table Hi-Table Hill give the input SNR thresholds of the beamformers in the three scenarios. From (l22l ). 
the input SNR thresholds can be determined as the following equation 

N 71 



SNR 



TO 



(48) 



L N — + 71) 

The values of the thresholds given in the tables are in accord with what are shown in the corresponding 
figures in the same scenarios. The thresholds of the proposed MIC-MPB scheme are far more less than 
those of the Maximin or PAPC scheme, and remain constants when INRs increase. The thresholds of the 
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TABLE III 

Input SNR Thresholds of the beamformers in five tones case 



Matrix Pair 


Input SNR Thresholds SNR T0 (dB) 


Beamformers 


INR = 10 dB 


INR = 20 dB 


INR = 30 dB 


MIC-MPB 


-0.64 


-0.56 


-0.55 


Maximin 


16.4 


26.4 


36.4 


PAPC 


oo 


oo 


oo 



20 
10 

-10 
-20 

m -30 

T3 

S - 50 



-90 
-100 
-110 

" 12 -?0 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90 

Direction of Arrival (degree) 

Fig. 6. The array patterns corresponding to the MIC-MPB, Maximin and PAPC with SNR = 10.9 dB and INR = 30 dB in 
two periodically repeated white noise case. 

Maximin beamformer increase the same amount accordingly when INRs increase 10 dB. The thresholds 
of the PAPC beamformer also show its failure because the values are always infinity in the three scenarios. 

Fig. [6] and Fig. [7] demonstrate the array patterns of the MIC-MPB, Maximin, and PAPC beamformer in 
the two periodically repeated white noise case. In Fig. |6j the proposed MIC-MPB scheme can correctly 
receive the SOI and null the interferes, but the Maximin or PAPC beamformer receives the interferers 
and forms a side-lobe in the direction of the SOI. The figure indicates that the MIC-MPB works at the 
operating area while both the Maximin and PAPC beamformer work at the failure area for SNR = 10.9 
dB and INR = 30 dB (c.f. Fig. [3] and Table In Fig. |7J the received signal power is very large and 
SNR = 40.9 dB is much larger than SNRjo of the MIC-MPB and Maximin algorithm, so both algorithms 



Interferer 1 the SOI Interferer 2 




Maximin 

— MIC-MPB 

— PAPC 



October 4, 2010 



DRAFT 



20 



20 



Interferer 1 



the SOI 



Interferer 2 



Maximin 

— MIC-MPB 

— PAPC 




10 




-70 



-80 



-90 



10 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90 

Direction of Arrival (degree) 



Fig. 7. The array patterns corresponding to the MIC-MPB, Maximin and PAPC with SNR = 40.9 dB and INR = 30 dB in 
two periodically repeated white noise case. 



Fig. 8. Geometrical interpretation of different thresholds of the MIC-MPB, Maximin and PAPC beamformer. 

can work properly. However, the Maximin beamformer just form a side-lobe or a shallow notch in the 
direction of the interferers. For the PAPC beamformer, a very deep null are placed in the direction of 
the SOI for (3 / and Rj contains part of the SOI, which can partly explain why G decreases when 
SNR increases shown in the above figures. 

Geometrical interpretation of different thresholds of the MIC-MPB, Maximin and PAPC beamformer 
in the scenarios can be illustrated by Fig. [8] For the beamformer with one interference channel or one 
dimensional interference subspace X, the condition (l2Tb which make 71 bounded is equivalent to the 
condition that requires the projected vectors of Hj and h.5 onto V/ must be in one line (c.f. Fig. [8). 
But this condition can hardly be satisfied for uncertainty of the characteristics of the interferers. For the 
proposed beamforming scheme with multiple interference channels, since there are multiple base vectors 
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Two Interference 
Channels 



V, = 7?{s,} is the subspace 
spanned by the waveforms of the interferers in one period 
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Fig. 9. Normalized average output SINR corresponding to the MIC-MPB, Maximin, and PAPC-RLS vs. sample size under 
various SNR in 1000 trials. 

in the interference channel, the condition can be easily satisfied. 

B. Performance of convergence rate with finite samples 

In this subsection, we compare the performance of convergence rate of the MPBs with finite samples. 
In the simulations, we assume the receiver has an array of ten elements (L = 10) with half wavelength 
spacing, and receives a single path SOI from 20°. There are seven MAIs, with INR of 40 dB and DOAs 
of 35°, -35°, -45°, 0°, -50°, -60° and 45°, respectively. Moreover, a broadband BPSK jamming also 
arrives from 60° with INR of 40 dB. These parameters have been verified not to cause obvious threshold 
effects of the Maximin and PAPC beamformer. Since there are two different approaches-stochastic gradient 
method O and recursive least squares (RLS) method ifTOll for PAPC beamformer to search the optimal 
weight vector in the literature, we name the algorithms as PAPC-SG and PAPC-RLS respectively for 
notational convenience. Fig.[9]shows the normalized output SINRs, defined as the ratio of output SINRs to 
the optimum value SINR opt under given SNR, which are calculated by averaging over 1000 independent 
trials. We observe that the proposed MIC-MPB scheme converges to the optimum performance within a 
few symbols, and is independent of the desired signal strength. In contrast, the PAPC-RLS and Maximin 
schemes require much more symbols and the performance of PAPC-RLS degrades when the input SNR 
increases. These results confirm the performance improvement of the MIC-MPB scheme, which extracts 
more effective samples per data symbol and eliminates the desired component in interference subspace 
1. 
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Fig. 10. Output SINR vs. time of the four beamformers in dynamic multiple access channel in 1000 trials. 



We also simulate the performance of different adaptive algorithms for dynamic multiple access channels. 
In this simulation, the input SNR is fixed to 20 dB. Fig. [10] compares the tracking ability of the PAPC- 
SG, PAPC-RLS, Maximin algorithm and MIC-MPB algorithm presented in section [TV} The DOAs of 
the seven MAIs are identical to the previous simulation. The first two MAIs are 8 dB stronger than the 
power of the SOI and the others are 40 dB stronger. The time they enter the channel are marked in the 
figure. The results demonstrate that the proposed recursive algorithm can null the new interferers within 
a few symbols, much faster than the other three algorithms. 

C. Performance when there are multipaths with identical delays 

In practice, the scatterers local to the mobile will cause an angular spread of about 3° at a distance 
of 1 km [22], and the relative delays between the multipaths are generally small. Thus, the assumption 
that the relative delays are greater than one chip may not hold. In this subsection, we will show that the 
proposed beamformer still work well under such condition. 

Assume there are D{ paths for the ith user. We first define a set Ui = {1, 2, . . . , D{} = (J s =i so 
that the subset Wj iS satisfies 

1) Vs ^ a', u i;S n Ui, s > = 0; 

2) Vj, j £ Ui tS , riij = riiji = n.j s . 

where riij, riiy, and all denote the equivalent propagation delays of certain paths. Thus, Ui jS contains 
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Fig. 11. Array patterns of MIC-MPB for paths with different delays. 



all the zth user's path indices of the same delay. As a result, we can rewrite (@]) as 

A/-1 S ( +oo Q 

x ( n ) =J2J2Y1 bi(k)ci(n - n is - kN)a(9 ls ) + ^ z q (n)a(6 q ) + v(n) (49) 

i=0 s=l k=—oo q=l 

where a(f?j s ) = Ylj^u ■y^ij a (^ij) * s tri e compound steering vector. For the desired user (i = 0) and 
Vj G Uo,s, the matrices and Rj will only depend on s, so we denote them as R^ and Rr iS 
respectively. The sth beamformer is then 

w opt , s = nR^^a.(e 0s ) =n y2 V^i R x, 1 s a ( 6 'oi)' ( 50 ) 

which means that the sth beamformer will cancel all other signals except the ones having the delay of 
uq s . Moreover, multiple beams will be formed to collect and combine the multipath components from 
different directions. Therefore, the algorithm is still applicable in such situation, and the only variation 
is that just So beamformers are required. 

Fig. [TT] and Fig. [12] show the simulated array patterns when the delays are different (dash line and dot 
lines), and the array pattern when the delays are identical (solid line). In the simulation, array elements 
L = 10 with half wavelength spacing are considered. Two users (M = 2) communicates with the receiver. 
The first user is the desired one and the second user acts as an MAI. There is a BPSK jammer from 
40° and INR = 40 dB. The bandwidth of the broadband jammer is 1/T C . Each user has two paths with 
equal power. The DOAs of the two desired paths are 0° and 12°. The paths of the second user arrive 
from —10° and —50°, and are 20 dB stronger than each path of the desired user. The input SNR for 
each desired path is 15 dB. In the former situation, the proposed MIC-MPB scheme forms two different 
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Fig. 12. Array patterns of MIC-MPB for paths with different and identical delays. 

beams to collect the two paths respectively, and each beamformer will suppress the other path besides the 
MAIs and the jammer. If the two desired paths have the identical delay, then one uniform beam will be 
formed to receive them, only nulling the MAIs and the jammer. Fig. [12] also shows when delays are not 
discriminable within one-chip period, two different beams will still be formed, but the two desired path 
are both collected by each beam. This implies that the proposed approach is robust to angular spread, 
where the delay spread is small. 

VI. Conclusion 

In this paper, we presented the principles for designing the projection space which are closely correlated 
with the ability of suppressing structured interference and system finite sample performance. According 
to the principles, we proposed an MIC-MPB scheme for CDMA systems which can be efficiently 
implemented by FFT. We also derived an adaptive algorithm for the beamformer. Computation and 
simulation results show that the proposed beamformer has a small and bounded SNR threshold, and 
can achieve the optimum SINR regardless of the received power of interference in the scenarios with 
structured interference. Furthermore, the various simulation results illustrate that the proposed MIC-MPB 
scheme has better finite sample performance, faster convergence rate and more superior tracking capability 
in the dynamical environment than the existing MPBs. 
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